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Abstract 

We examine the transition of a particle across the singularity of the compactified Milne (CM) 
space. Quantization of the phase space of a particle and testing the quantum stability of its 
dynamics are consistent to one another. One type of transition of a quantum particle is described 
by a quantum state that is continuous at the singularity. It indicates the existence of a deterministic 
link between the propagation of a particle before and after crossing the singularity. Regularization 
of the CM space leads to the dynamics similar to the dynamics in the de Sitter space. The CM 
space is a promising model to describe the cosmological singularity deserving further investigation 
by making use of strings and membranes. 
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I. INTRODUCTION 



Presently available cosmological data indicate that known forms of energy and matter 
comprise only 4% of the makeup of the Universe. The remaining 96% is unknown, called 
'dark', but its existence is needed to explain the evolution of the Universe 0, The dark 
matter, DM, contributes 22% of the mean density. It is introduced to explain the observed 
dynamics of galaxies and clusters of galaxies. The dark energy, DE, comprises 74% of the 
density and is responsible for the observed accelerating expansion. These data mean that 
we know almost nothing about the dominant components of the Universe! 

Understanding the nature and the abundance of the DE and DM within the standard 
model of cosmology has difficulties 0, 1^ . These difficulties have led many physicists to seek 
anthropic explanations which, unfortunately, have little predictive power. An alternative 
model has been proposed by Steinhardt and Turok (ST) 0, 0, 01- It is based on the idea 
of a cyclic evolution, CE, of the Universe. The ST model has been inspired by string/M 
theories 0]. In its simplest version it assumes that the spacetime can be modelled by the 
higher dimensional compactified Milne, CM, space. The attraction of the ST model is that 
it potentially provides a complete scenario of the evolution of the universe, one in which the 
DE and DM play a key role in both the past and the future. The ST model requires DE for 
its consistency, whereas in the standard model, DE is introduced in a totally ad hoc manner. 
Demerits of the ST model are extensively discussed in |9i]. Response to the criticisms of 
can be found in 

The mathematical structure and self-consistency of the ST model has yet not been fully 
tested and understood. Such task presents a serious mathematical challenge. It is the 
subject of our research programme. 

The CE model has in each of its cycles a quantum phase including the cosmological 
singularity, CS. The CS plays key role because it joins each two consecutive classical phases. 
Understanding the nature of the CS has primary importance for the CE model. Each CS 
consists of contraction and expansion phases. A physically correct model of the CS, within 
the framework of string/M theory, should be able to describe propagation of a p-brane, 
i.e. an elementary object like a particle, string and membrane, from the pre-singularity to 
post-singularity epoch. This is the most elementary, and fundamental, criterion that should 
be satisfied. It presents a new criterion for testing the CE model. Hitherto, most research 
has focussed on the evolution of scalar perturbations through the CS. 

Successful quantization of the dynamics of p-brane will mean that the CM space is a 
promising candidate to model the evolution of the Universe at the cosmological singularity. 
Thus, it could be further used in advanced numerical calculations to explain the data of 
observational cosmology. Failure in quantization may mean that the CS should be modelled 
by a spacetime more sophisticated than the CM space. 

Preliminary insight into the problem has already been achieved by studying classical and 
quantum dynamics of a te st p article in the two-dimensional CM space The present 

paper is a continuation of [10| and it addresses the two issues: the Cauchy problem at the 
CS and the stability problem in the propagation of a particle across the CS. Both issues 
concern the nature of the CS. 

In Sec. II we define and make comparison of the two models of the universe: the CM 
space and the regularized CM space. The classical dynamics of a particle in both spaces is 
presented in Sec. III. The quantization of the phase space of a particle is carried out in Sec. 
IV. In Sec. V we examine the stability problem of particle's dynamics both at classical and 
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quantum levels. We summarize our results, conclude and suggest next steps in Sec. VI. 



II. SPACETIMES 
A. The CM space 

For completeness, we recall the definition of the CM space used in 0]. It can be specified 
by the following isometric embedding of the 2d CM space into the 3d Minkowski space 



y^{t,e) =tVl + r2, y\t,e) =rts\Ti{e/r), y\t , 6) = rt cos{e / r) , (1) 
where {t, 0) G x §^ and < r G is a constant labelling compactifications . One has 



{yr-iyr-iy'r = o. (2) 



1 + 



Eq. presents two cones with a common vertex at {y'^,y^,y'^) = (0,0,0). The induced 
metric on ^ reads 

ds^ = -de + edd\ (3) 

Generalization of the 2d CM space to the Nd spacetime has the form 

ds^ = -df + dx^dxk + t^de"^, (4) 

where t, x'^ G M^, ^ G (fc = 1, . . . , — 2). One term in the metric (jU disappears/appears 
at t = 0, thus the CM space may be used to model the big-crunch/big-bang type singularity 
In what follows we restrict our considerations to the 2d CM space. Later, we make 
comments concerning generalizations. 

It is clear that the CM space is locally isometric to the Minkowski space at each point 
except the vertex t = 0. The CM space is not a manifold, but an orbifold due to this vertex. 
The Riemann tensor components vanish for t 7^ and cannot be defined at t = 0, since one 
dimension disappears/appears there. There is a space-like singularity at t = of removable 
type because any time-like geodesic with t < can be extended to some time-like geodesic 
with t > 0. However, such an extension cannot be unique due to the Cauchy problem for the 
geodesic equation at the vertex (compact dimension shrinks away and reappears at t = 0). 



B. The RCM space 

Since trajectory of a test particle coincides (by definition) with time-hke geodesic, there is 
no obstacle for the test particle to reach and leave the CS. However, the Cauchy problem for 
a geodesic equation at the CS is not well defined. As the result, a test particle 'does not know 
where to go' at the singularity. Thus, the singularity acts as 'generator' of uncertainty in the 
propagation of a test particle from the pre- singularity to post-singularity era. In the present 
paper we propose to solve this problem by replacement of a test particle by a physical one. 
The test and physical particles differ in a number of ways. For instance, physical particle's 
own gravitational field effects its motion fll| and may modify the singularity of the CM 
space. We assume that these effects may be modelled by replacing the CM space by a 
regularized compactified Milne, RCM, space in such a way that the big-crunch/big-bang 
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type singularity of the CM space is replaced by the big-bounce type singularity. In the 
RCM space the Cauchy problem does not occur because compact space dimension does not 
contract to a point, but to some 'small' value. As the result the propagation of a particle is 
uniquely defined in the entire spacetime. Particle's prop agation in the RCM space is similar 
to the corresponding one in the de Sitter space l^, 13| . 

We define the RCM space by the following embedding into the 3d Minkowski space 



and we have the relation 



rVW+l^sm{e/r), y^{t,e) = cos{e/r), (5) 



(y'r - {y'f - if) 



,2\2 



1 + r2 

The induced metric on the RCM space reads 



2 2 

— e r . 



dsj = -(1 + 



t' + e' 



(6) 



(7) 



where e G M is a small number. It is clear that now the space dimension 9 does not shrink 
to zero at i = 0. The scalar curvature has the form 



7^. 



(e2(l + r2) +t2)2 



(8) 



and the Einstein tensor corresponding to the metric ((7j) is zero, thus ((Tj) defines some vacuum 
solution to the 2d Einstein equation. 
It is evident that at t 7^ we have 



\im.dsl = -dt^ + t^de'^ 



and 



lim 71^ 



0. 



It is obvious that © turns into (0) as e — > 0. 



(9) 





Figure 1: Embeddings of RCM and CM spaces. 



Figure 1 presents the RCM and CM spaces embedded into the 3d Minkowski space. We 
can see that the big-crunch/big-bang singularity of the CM space is represented in the RCM 
by the big-bounce type singularity. 
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III. CLASSICAL DYNAMICS 



An action integral, A, describing a relativistic test particle of mass m in a gravitational 
field gki, {k, / = 0, 1) may be defined by 

/• k • I 
drL{r), L{T):=-{^g,,-e), (10) 

where x'^ := dx^/dr, t is an evolution parameter, e(r) denotes the 'einbein' on the world- 
line, x° and x^ are time and space coordinates, respectively. 
In case of the CM and RCM spaces the Lagrangian reads 

L.ir) = ^Jie + e^)e^-il + ^)i^-e^), (11) 

where e = corresponds to the CM space. The action (fTUI) is invariant under reparametriza- 
tion with respect to r. This gauge symmetry leads to the constraint 

:= PkPi g'^ + m' = - -A^ + = 0, (12) 

\t ^ e ) 1 + j2_,_^2 

where := dL^/dt and pe '■= dL^/dO are canonical momenta, and where g'^^ denotes an 
inverse of the metric gui defined by the line element (jZ)) (case e = corresponds to the CM 
space). 

Variational principle applied to (fTUI) gives equations of motion of a particle 

d ^ d dL ^ dL ^ 

Since during evolution of the system po is conserved, due to (fT^ . we can analyze the 
behaviour of pt by making use of the constraint (fT^ . In case of the CM space (e = 0), for 

7^ there must be pt oo as t ^ 0. This problem cannot be avoided by different choice 
of coordinates^. It is connected with the vanishing/ appearance of the space dimension 9 at 
t = 0. Another interpretation of this problem is that different geodesies cross each other 
with the relative speed reaching the speed of light as they approach the singularity at t = 0. 

The dynamics of a physical particle in the RCM space (e 7^ 0) does not suffer from such a 
problem, since for 7^ the momentum component pt does not need to 'blow up' to satisfy 

A. Geodesies in CM and RCM spaces 

It was found in ^] an analytic general solution to (jl3|) . for e = 0, in the form 

^^(t) = ^^o-sinh-(£), (14) 



^ The system of coordinates we use, (i, 6') e x is natural for the spacetimes with the topologies 
presented in Fig. 1. 
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where {pe, 9o) eM} xE>^. It is clear that geodesies (fT^ 'blow up' at t = 0, which is visualized 
in Fig. 2. 

For e 7^ 0, Eqs. (0 read 

!I^^l±^ =p^ = const, e' = {l + - if + (15) 

e + 

and 

(i + i^^)*--(i + i^^)(;)*-^' + ''=*-o. (16) 



From (jl5|) and we get 



.dtJ w?{t^ + + plif^ + e^y 
where G M^. General solution to (fT7|) reads 



(17) 



e[t) = »o + P« f ^ dr^ ro'lj^^^^^^ (18) 

where Oq G The integral in (jlSp cannot be calculated analytically. Numerical solution of 
(fT7|) is presented in Fig. 2. 




Figure 2: Geodesies in ROM space (the left graph) and CM space (the right graph). 



Figure 2 shows that a geodesic in the RCM space is bounded and continuous in the 
neighborhood of the singularity. In contrary, a geodesic in the CM space (drawing by 
making use of (jl4|) ) blows up as t ^ ±0. 



B. Phase space and basic observables 

We define a phase space to be the set of independent parameters (variables) defining all 
particle geodesies. Thus the pase space, F, for (fTH|) reads 

r := {((T,p<,) I CT G mod 27rr, p„ G M}} = x (19) 

The Cauchy problem at the singularity results from the vanishing/appearance of the space 
dimension at t = 0. It is fairly probable that any simple regularization of the singularity 
of the CM space that prevents such collapse will lead to the cylindrical phase space (jl9j) . 
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In we have analyzed four types of propagations of a particle in the CM space. Now 
we can see that the regularization prefers the propagation in the CM space of the de Sitter 
type (see, Sec. Ill D of |10|]), because only in this case the phase space topology has the 
form dTTHl . 

Now, let us identify the basic canonical functions on the phase space, i.e. observables 
that can be used to define any composite observable of the underlying classical system. In 
case a phase space includes a variable with non-trivial topology, i.e. different from M^, it is 
a serious problem. However, it has been solved in two (equivalent) ways not long ago. In 
what follows we use the method used in the group theoretical quantization (see, and 
references therein). In the next section we explain relation with another method. 

A natural choice of the basic functions on ()19|) is 

S* := sin(cr/r), C := cos(cr/r), P := rp„. (20) 

The basic observables 5* and C are smooth single- valued functions on (contrary to cr). 
The observables satisfy the Euclidean algebra e(2) on T 

{S,C} = 0, {P,S} = C, {P,C} = -S, (21) 

where 

d ■ d • d • d • 

^ ' ^ ' dp^ da da dp„ ' ^^^"^ 



It is shown in that the Euclidean group E{2) can be used as the canonical group of 
the phase space T. 



IV. QUANTIZATION OF PHASE SPACE 

By quantization we mean finding an irreducible unitary representation of the symmetry 
group of the phase space of the underlying classical system. 

The group E{2) has the following irreducible unitary representation [l^ 

[U{a)ilj]{(3) := ^p[{f3 - a) mod 2ti], for rotations z e^"z, (23) 

[U{t)ip]{P) := [exp —i{a cos (3 + b sin for translations z^z + t, (24) 

where z = \z\ e*^, t = a + bi, and where ip G L'^{E>^). 

Making use of the Stone theorem, we can find an (essentially) self-adjoint representation 
of the algebra (j^H). One has 

[C,S]=0, [P,S] = -tC, [P,C]=tS, (25) 

where 

P^iP) := -^^^iP), S^iP) ■■= C^iP) := cos/3 ^(P). (26) 

The domain, Qx, of operators P, S, C reads 

■= e L\n') I ^ e C~[0, 2n], (^(")(27r) = e^V^^Ho), n = 0, 1, 2, . . . }, (27) 
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where < A < 27r labels various representations of e(2) algebra. The space Q\ is dense in 
L^(S^) so the unbounded operator P is well defined. As the operators S and C are bounded 
on the entire L^(S^), the space fix is a common invariant domain for all operators and their 
products. 

In [l^ we have found that the representation of the algebra specific to the case considered 
there in Sec. Ill D, has the form 

[a,U]=U, (28) 

with 

:= -^^m, MP) ■■= eMP), (29) 

where < /3 < 27r and (p G Qx. However, both representations, and (|^. are in fact 
the same owing to 

e'^ = cos /? + z sin p, [cos /3, sin p] = 0. (30) 

The space fix, where < A < 2tt, may be spanned by the set of orthonormal eigenfunc- 
tions of the operator a 

fk,x{P)-=i'^^y'^^eWiPik + \/27i), k = 0,±l,±2,... (31) 

However, the functions ()31|) are continuous on only in the case when A = 0. Thus, the 
requirement of the continuity removes the ambiguity of quantization. 



V. STABILITY OF SYSTEM 

To examine the stability problem of our system, we use the Hamiltonian formulation of 
the dynamics of a particle. 



A. Classical level 

By stability of the dynamics of a classical particle we mean such an evolution of a particle 
that can be described by the canonical variables which are bounded and continuous functions. 

Direct application of the results of jl^l gives the following expression for the extended 
Hamiltonian jiol. Il7|. H^, of a particle 

He = la (32) 

where is an arbitrary function of an evolution parameter r, and where $e is the first-class 
constraint defined by (fT^ . The equations of motion for canonical variables {t,6;pt,pe) read 



^ = {t,H.} = -air) Pf^l'^ , (33) 
dr t'^ + + r^e^ 



|.{«,^.}.C.M^, (34) 
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^ = {p,,i/,} = 0, (36) 



where 

, d-d- d- d- d- d- d- d- 

dt dpt dpt dt 86 dpe dpe dO 
To solve we use the gauge t = t. In this gauge (jH^ leads to 

j-2 _|_ ,2 _|_ »,2,2 

CM - -^J^- (37) 

Insertion of into and taking into account (fT^ gives 

d_ 2 ^ 2tpye'' 2tp^(t^ + 6^ + rV) 

ci/* (t2 + e2)(^2 + g2 + ^2e2) (t2 + e2)3 " ^ > 



2 _ / cipI \ + + r^e^ 



Solution to (pH|l reads 

where pe does not depend on time due to This result is consistent with the constraint 
flT^ if we put ci = 1 and C2 = rn? . Insertion of (jHIIj) into (jTTj) yields an explicit expression 
for Ce(t). Next, insertion of so obtained Ceit) into (jMj) gives (fTTj) with the solution (fTHjl . 
Thus, we have found complete solutions^ to the equations (jH^ - flHUj) . 

It results from the functional form of solutions that for e 7^ the propagation of canonical 
variables is regular, i.e. has no singularities for any value of time. One may also verify that 
the constraint equation (fT^ . with pt determined by (j^ . is satisfied for each value of time 
either. Since Ce(t), determined by P7jl is bounded, the Hamiltonian defined by (jH^ is 
weakly zero due to the constraint (fT^ . 

Now, let us analyze the case e = 0, which corresponds to the evolution of a particle in 
the CM space. The solutions of in the gauge t = t, are the following 

Co{t) = -1/pt, (40) 

p't=pl/t' + m' (41) 

9(t) = -sinh~^{pe/mt^ + const, (42) 

(where pe = const), in agreement with (fT^ . The constraint equation (fT2j) reads 

%=pl/t^-p',+m^ = 0. (43) 

It results from ()4H) and ()42|) that only for pg = the propagation of a particle in the CM 
space is regular for any value of time. The Hamiltonian (jH^ is also regular and is weakly 
equal to zero. Quite different situation occurs in the case pg 0. The equations (PT|) - (P3j) 
and the Hamiltonian are singular at t = 0. Thus the dynamics of a particle is unstable^. 



^ This way we have also verified, in the gauge t — t, the equivalence between our Hamiltonian and the 

Lagrangian formulations of the dynamics of a particle. 
^ The division of the set of geodesies into regular and unstable depends on the choice of coordinates, but it 

always includes the unstable ones. 
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B. Quantum level 

By stability of dynamics of a quantum particle we mean the boundedness from below of 
its quantum Hamiltonian. 

To construct the quantum Hamiltonian of a particle we use the following mapping (see, 
e.g. [13) 

pmg'^ □ := (-g)-y^d,[i-gy/'g'^di], (44) 

where g := det[gki] and dk '■= d/dx^. The Laplace-Beltrami operator, □, is invariant 
under the change of spacetime coordinates and it leads to Hamiltonians that give results 
consistent with experiments ^3 5 which has been used in theoretical cosmology (see, 
(3^ and references therein). 

In the case of the CM space the Hamiltonian, for t < or t > 0, reads 

The operator H was obtained by making use of (jUj), Q and (jH^ in the gauge = 2 (for 
e = 0). In this gauge^ the Hamiltonian equals the first class constraint Thus the Dirac 
quantization scheme (3, [l^ leads to the equation 

H^{e,t) = 0. (46) 

The space of solutions to pUj) defines the domain of boundedness of H from below (and 
from above). 

Let us find the non-zero solutions of (jlUj) . Separating the variables 



leads to the equations 
and 



^lj{e,t) := A{e) Bit) (47) 

d'^A/de'^ + p^A = 0, p G M (48) 
d'^B 1 dB mH'^ + p2 

where p is a constant of separation. Two independent continuous solutions on §^ read 

Ai{p,6) = aiCos{p9), A2{p,0) = a2sm{p6), 01,02 G M. (50) 

Two independent solutions on R (for t < or t > 0) have the form 

B^{p, t) = hi^J{ip, mt), B2{p, t) = b2^Y{ip, mt), 61, &2 e C, (51) 

where 3? J and are the real parts of Bessel's and Neumann's functions, respectively. Since 
p G M, the number of independent solutions is: 2 x 2 x 00 ( for t < and t > 0). 



^ In the preceding subsection, concerning the classical dynamics, the choice of gauge was different. But 
since the theory we use is gauge invariant, the different choice of the gauge does not effect physical results. 
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At this stage we define the scalar product on the space of solutions (fSUj) and (fFIj) as 
follows 

< 4'i\4'2 >■= dfJ' i^i 4'2, dfi := dd dt = \t\ dO dt, (52) 

where F := [— T, 0[x§^ (with T > 0) in the pre-singulaity epoch, and F :=]0,T] x §^ in the 
post-singularity epoch. We assume that the CM space can be used to model the universe 
only during its quantum phase, which lasts the period [— T, T] . 

Now we construct an orthonormal basis, in the left neighborhood of the cosmological 
singularity, out of the solutions (jSUj) and (j^D). One can verify that the solutions (jSUj) are 
orthonormal and continuous on §^ if ai = tt^^/^ = 02 and rp = 0, ±1, ±2, . . .. (This set 
of functions coincides with the basis (j31|) that spans the subspace Qx if we replace k by rp.) 
Some effort is needed to construct the set of orthonormal functions out of ^J{ip,mt) and 
^Y{ip,mt). First, one may verify that these functions are square-integrable on the interval 
[— T, T]. This is due to the choice of the measure in the scalar product which leads 

to the boundedness of the corresponding integrants. Second, having normalizable set of 
four independent functions, for each p, we can turn it into an orthonormal set by making 
use of the Gram-Schmidt procedure (see, e.g. jl^]). Our orthonormal and countable set of 
functions may be used to define the span J^. The completion of JF in the norm induced by 
the scalar product defines the Hilbert spaces L^(F x S^,dp). It is clear that the same 
procedure applies to the right neighborhood of the singularity. 

p(e,t) p(e,t) 




(a) near the singularity (b) on large scale 

Figure 3: Probability density corresponding to '0(6', = Ai{0,9) 5RJ(0,t) 

Finally, we can prove that the Hamiltonian ()45|) is self-adjoint on L^(F x S>^,dp). The 
proof is immediate if we rewrite ()46|) in the form 

□ 7/; = -m^ ij. (53) 

It is evident that on the orthonormal basis that we have constructed above the operator □ 
is an identity operator multiplied by a real constant — m^. The operator □ is bounded since 

||n|| := sup !!□ ipW = sup II — ipW = w? < 00, (54) 

m=i 11^11=1 

where jl^/'jl := ^J < > . The operator □ is also symmetric, because m is a real constant. 
Since □ is bounded and symmetric, it is a self-adjoint operator (see, e.g. jl^). Clearly, the 
self-adjointness of the Hamiltonian ()45p results from the self-adjointness of □. 
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p(e,t) P(e,t) 




(a) near the singularity (b) on large scale 

Figure 4: Probability density corresponding to tl^{9,t) = ^i(0,^) ^Y(0,t) 

We have constructed the two Hilbert spaces: one for the pre- singularity epoch, 'H^~\ and 
another one to describe the post-singularity epoch, Ti.^^^ Next problem is to 'glue' them 
into a single Hilbert space, TC = L^([— T, T] x S^,dfi), that is needed to describe the entire 
quantum phase. From the mathematical point of view the gluing seems to be problematic 
because the Cauchy problem for the equation (pUj) is not well defined^ at t = 0, and because 
we have assumed that t 7^ in the process of separation of variables to get Eqs. (jlSjl and 
fj49|) ■ However, arguing based on the physics of the problem enables the gluing. First of all 
we have already agreed that a classical test particle is able to go across the singularity (see, 
subsection H B). One can also verify that the probability density 

P{e,t) := = \t\ mO,t)\' (55) 

is bounded and continuous in the domain [— T, T] x Figures 3 and 4 illustrate the 
behavior of P{6,t) for two examples of gluing the solutions having p = 0. The cases with 
p ^ have similar properties. Thus, the assumption that the gluing is possible is justified. 
However one can glue the two Hilbert spaces in more than one way, as it was done in the 
quantization of the phase space in our previous paper In what follows we present two 
cases, which are radically different. 

1. Deterministic propagation 

Among all solutions (|51|) there is one, corresponding to p = 0, that attracts an attention 
It reads 

Bi{0, mt) = hi 3f?J(0, mt), hi E M, (56) 
and has the following power series expansion close to t = 

2 4 fi 

Bi{<d,x)/hi = l- — + - — + 0\x% (57) 

IV , ;/ 1 4 64 2304 ^ ^ ^ ' 

It is visualized in Fig. 5a. The solution (j^UI) is smooth at the singularity, in spite of the fact 
that (jl^ is singular at t = 0. 



^ Except one case discussed later. 
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Figure 5: Continuous (a) and singular (b) propagations of a particle with p = 0. 

It defines a solution to (pUj) that does not depend on 9, since the non-zero solution (p])) 
with p = is just a constant. Thus, it is unsensitive to the problem that one cannot choose 
a common coordinate system for both t < and t > 0. 

The solution Bi (and the trivial solution Bq := 0) can be used to construct a one- 
dimensional Hilbert space H = L^{[—T,T] x E>^,dp). The scalar product is defined by (fH^ 
with r replaced by F := [— T, T] x It is obvious that the Hamiltonian is self-adjoint on 

n. 

The solution (jSUj) is continuous at the singularity. It describes an unambiguous propaga- 
tion of a quantum particle. Thus, we call it a deterministic propagation. It is similar to the 
propagation of a particle in the RCM space considered in the next subsection. 

Since is a second order differential equation, it should have two independent solutions. 
However, the second solution cannot be continuous at t = 0. One may argue as follows: The 
solution (j56|) may be obtained by ignoring the restriction t ^ and solving (j49|) with the 
following initial conditions 

5(0,0) = 1, dB{0,0)/dt = 0. (58) 

Equations (j49|) and (j58j) are consistent, because the middle term of the r.h.s. of (j49|) may be 
equal to zero due to so the resulting equation would be non-singular at t = 0. Another 
initial condition of the form 5(0,0) = const and dB{0,0)/dt = would be linearly 
dependent on (j3Hj) . Thus, it could not lead to the solution which would be continuous at 
t = and linearly independent on ()56|) . 

This qualitative reasoning can be replaced by a rigorous derivation using the power series 
expansion method Applying this method one obtains that near the singularity t = 
the solution to ()49p behaves like and that the corresponding indicial equation reads 

= -p\ (59) 

Thus, for p 7^ the two solutions behave like t^*'' , i.e. are bounded but not continuous 
(see, Eq. (jSH))- For p = the indicial equation has only one solution = which leads 
to an analytic solution to (PU)) defined by (jSHI)- In such a case, it results from the method 
of solving the singular linear second order equations ^], the second solution to (j49j) may 
behave like In |t| . In fact it reads jl^ 

^2(0, mt) = 62 3ftF(0, mt), 62 G M, (60) 
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and is visualized in Fig. 5b. It cannot be called a deterministic propagation due to the 
discontinuity at the singularity t = 0. 



2. Indeterministic propagation 

All solutions (jSU, except (fSUI) . are discontinuous at t = 0. This property is connected 
with the singularity of at t = 0. It is clear that due to such an obstacle the identification 
of corresponding solutions on both sides of the singularity is impossible. However there are 
two natural constructions of a Hilbert space out of 'H^~'^ and H^^^ which one can apply: 
(a) Tensor product of Hilbert spaces 

The Hilbert space is defined in a standard way as H := n^-^ ® and it consists of 
functions of the form 

fih,9,;t2,e,) ^ (/(-) ® /W)(ti,^^i;t2,^2) := f^-\h,e,) f^+\t^,e^), (ei) 

where /'•^•* G TY*-^-* and /^+-' G TY'-^^ The scalar product reads 

< / I ^ >:=< /^-^l 9^-^ > < f^^-^l 9^^^ >, (62) 



where 



and 



< 9^-^ >■■= f dt^ r dOr /(-)(ti,^i) 9'^-\h,e,) (63) 

J-T Jo 



< /^+^| 9^^^ >■■= / dt2 / de^ \t2\ /(+Ht2,^2) 9^^\t2,e2). (64) 
Jo Jo 

The action of the Hamiltonian is defined by 

(/(-) ® /W) := (^/(-)) ® /(+) + /(-) ® (jy/ (+)) . (65) 

The Hamiltonian is clearly self-adjoint on Ti. 

The quantum system described in this way appears to consist of two independent parts. 
In fact it describes the same quantum particle but in two subsequent time intervals separated 
by the singularity at t = 0. 
(b) Direct sum of Hilbert spaces 

Another standard way \2l\ of defining the Hilbert space is Ti. := IHS^^ @'H^^\ The scalar 
product reads 

< /1I/2 >:=< fV\ft^ > + < fi^^ft^ >, (66) 

where 

fk := ift\fl^^) G n^-^ X T^W, k = l,2, (67) 

and where fjr^ and fj.'^^ are two completely independent solutions in the pre-singularity and 
post-singularity epochs, respectively. (The r.h.s of is defined by and 
The Hamiltonian action on 7i reads 

n 3 (/(-), /W) ^ H{f^'\f^^^) := {Hf^-\Hf^^^) G H. (68) 

It is obvious that H is self adjoint on Ti. 



14 



By the construction, the space 7i^~^ 0'^*-''''' includes vectors hke (/^~-*,0) and (0, Z*^^-'), 
which give non- vanishing contribution to (but yield zero in case (j6!2j) ). The former 
state describes the annihilation of a particle at t = 0. The latter corresponds to the creation 
of a particle at the singularity. These type of states do not describe the propagation of a 
particle across the singularity. The annihilation/creation of a massive particle would change 
the background. Such events should be eliminated from our model because we consider a 
test particle which, by definition, cannot modify the background spacetime. Since and 
being vector spaces, must include the zero solutions, the Hilbert space 0'^^*^^-* 
cannot model the quantum phase of our system. 



C. Regularization 

In the RCM space the quantum Hamiltonian, H^, for any t G [— T, T], reads (we use the 

gauge Ce = 2) 

_ Vt^ + + eV^ &^ _ e + e^ 

Since the Hamiltonian is equal to the first-class constraint, the physical states are solutions 
to the equation 

if,^ = 0. (70) 
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Figure 6: Properties of the symmetric solution with e = 0.5, /O = 1 and r = 1: (a) the solution 
corresponding to Bs{l, 0) = 1 and dBs{l, 0)/dt = as the initial values for (f7^ . (b) the probability 
density P{9,t) corresponding to ^i(l,0) and Bs{l,t). 



As in the case of e = 0, the space of solutions to (fTHjl defines the domain of boundedness 
of i^e- Substitution ^(^,t) = A{e)B{t) into ^ yields 



d9^ 



+ = 0, pe 



and 



(t2 + e2)2 ^25 t{t^ + e''){t^ + + 2rh^) dB 



t2 + e2 + r2e2 dt^ 



+ 



(t2 + e2 + 2r2e2)2 



dt 



+ m\t' + + p^)B = 0. 



(71) 



(72) 
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B,(l,t) P(0,t) 




(a) (b) 

Figure 7: Properties of the anti-symmetric solution with e = 0.5, p = 1 and r = 1: (a) the solution 
corresponding to Ba{l, 0) = and dBa{l,0)/dt = 1 as the initial values for (f7^ . (b) the probability 
density P{0,t) corresponding to Ai(l,0) and Ba{l,t). 



The equation (fTTj) looks the same as in the non-regularized case f)lHj) so the two indepen- 
dent solutions on read 

A, (p, 6) = 7r-i/2 cos(pe) , A2 (p, e) = 7r-i/2 sin(p^) , (73) 

where rp = 0, ±1, ±2, . . . (orthogonality and continuity conditions). 

Equation (f72|) is non-singular in [— T, T] so for each p it has two independent solutions 
which are bounded and smooth in the entire interval. One may represent these solutions by 
a symmetric and an anti-symmetric functions. We do not try to find the analytic solutions 
to (j72|) . What we really need to know are general properties of them. In what follows we 
further analyze only numerical solutions to (j72p by making use of jl9| . 

Since in the regularized case the solutions are continuous in the entire interval [— T, T], 
the problem of gluing the solutions (the main problem in case e = 0) does not occur at all. 
Thus, the construction of the Hilbert space. He, by using the space of solutions to (f7T|) and 
fl7^ is straightforward. The construction of the basis in He may be done by analogy to 
the construction of the basis in H^~^ (described in the subsection B). The only difference is 
that now t G [— T, T] and instead of we use the solutions to (f?^. Let us denote them 
by Bi{p,mt), where i = s and i = a stand for symmetric and anti-symmetric solutions, 
respectively. Figures 6 and 7 present two examples of solutions to (f7^ for p = 1 and the 
corresponding probability densities. We can see that P{9,t) is a bounded and continuous 
function on [— T, T] x as in the case of e = (cp with Figs. 3 and 4). 

We define the scalar product as follows 

< ^^1^2 >:= j d^i ^1 dp := dO dt = Vt^ + + r'^e^ dO dt, (74) 

where F := [— T, T] x S^, and where an explicit form of dp is found by making use of 

It is evident that He is self-adjoint on He- The main difference between the deterministic 
case with e = and the present case e > is that in the former case the Hilbert space is 
one dimensional (p = 0), whereas in the latter case it is 2 x 2 x 00 dimensional (rp G Z). 
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VI. SUMMARY AND CONCLUSIONS 



The Cauchy problem at the cosmological singularity of the geodesic equations may be 
'resolved' by the regularization, which replaces the double conical vertex of the CM space 
by a space with the vertex of the big-bounce type, i.e. with non-vanishing space dimension 
at the singularity. We have presented a specific example of such regularization of the CM 
space. Both classical and quantum dynamics of a particle in the regularized CM space are 
deterministic and stable. We have examined these aspects of the dynamics at the phase 
space and Hamiltonian levels. The classical and quantum dynamics of a p article in the 
regularized CM space is similar to the dynamics in the de Sitter space [ij, We are 
conscious that our regularization of the singularity is rather ad hoc. Our arguing (presented 
at the beginning of Sec. IIB) that taking into account the interaction of a physical particle 
with the singularity may lead effectively to changing of the latter into a big-bounce type 
singularity should be replaced by analyzes. However, examination of this problem is beyond 
the scope of the present paper, but will be considered elsewhere. 

The classical dynamics in the CM space is unstable (apart from the one class of geodesies). 
However, the quantum dynamics is well defined. The Cauchy problem of the geodesies is 
not an obstacle to the quantization. The examination of the quantum stability has revealed 
surprising result that in one case a quantum particle propagates deterministically in the 
sense that it can be described by a quantum state that is continuous at the singularity. This 
case is very interesting as it says that there can exist deterministic link between the data of 
the pre-singularity and post-singularity epochs. All other states have discontinuity at the 
singularity of the CM space, but they can be used successfully to construct a Hilbert space. 
This way we have proved the stability of the dynamics of a quantum particle. 

At the quantum level the stability condition requiring the boundedness from below of 
the Hamiltonian operator means the imposition of the first-class constraint onto the space 
of quantum states to get the space of physical quantum states. The resulting equation 
depends on all spacetime coordinates. In the pre-singularity and post- singularity epochs 
the CM space is locally isometric to the Minkowski space [lO|. Owing to this isometry, the 
stability condition is in fact the Klein-Gordon, KG, equation. The space of solutions to 
the KG equation in these two epochs and the corresponding Hilbert space are fortunately 
non-trivial ones, otherwise our quantum theory of a particle would be empty. 

Quantization of the phase space carried out in Sec. IV (and in our previous paper [l( 
corresponds to some extent to the method of quantization in which one first solves constraints 
at the classical level and then quantize the resulting theory. Quantization that we call here 
examination of the stability at the quantum level, is effectively the method in which we 
impose the constraint, but at the quantum level. The results we have obtained within both 
methods of quantization are consistent. It means that the quantum theory of a particle in 
the compactified Milne space does exist. The CM space seems to model the cosmological 
singularity in a satisfactory way^. 

It turns out that the time-like geodesies of our CM space may have interpretation in terms 
of cosmological solutions of some sophisticated higher dimensional field theories [i^ . 
We have already discussed some aspects of this connection in our previous paper (see Sec. 
5 of 10]). Presently, we can say that in one case (see Sec. 4 of jl^ and Sec. V.B.I of 



^ Our result should be further confirmed by the examination of the dynamics of a particle in a higher 
dimensional CM space. 
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this paper) this analogy extends to the quantum leveh transition of a particle through the 
cosmological singularity in both models is mathematically well defined. In both cases the 
operator constraint is used to select quantum physical states. Elaboration of this analogy 
needs an extension of our results to the Misner space (that consists of the Milne and Rindler 
spaces) because it is the spacetime used in 0|. Another subtlety is connected with the fact 
that we carry out analysis in the compactified space, whereas the authors of use the 
covering space. 

There exists another model to describe the evolution of the universe based on string/M 
theory. It is called the pre-big-bang model 13]. However, the ST model is more self- 
consistent and complete. 

Other sophisticated model called loop quantum cosmology, LQC, is based on non- 
perturbative formulation of quantum gravity called the loop quantum gravity, LQG 0, 
It is claimed that the CS is resolved in this approach However, this issue seems to be 
still open due to the assumptions made in the process of truncatin g th e infinite number de- 
grees of freedom of the LQG to the finite number used in the LQC [sol, l3l|- This model has 
also problems in obtaining an unique semi-classical approximations |32], which are required 
to link the quantum phase with the nearby classical phase in the evolution of the Universe. 
For response to i32] we recommend [s^. 

Quantization of dynamics of extended objects in the CM space is our next step. There 
exist promising results on propagation of a string and membrane (s^. Issj]. However, these 
results concern extended objects in the low energy states called the zero-mode states and 
quantum evolution is approximated by a semi-classical model. Recently, we have quantized 
the dynamics of a string in the CM space rigorously 's^, but our results concern only the 
zero-mode state of a string. For drawing firm conclusions about the physics of the problem, 
one should also examine the non-zero modes. Work is in progress. 
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